The impact of hydrodynamic interactions on the preferential concentration of inertial 

particles in turbulence 
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We consider a dilute gas of inertial particles transported by the turbulent flow. Due to inertia the 
particles concentrate preferentially outside vortices. The pair-correlation function of the particles' 
concentration is known to obey at small separations a power-law with a negative exponent, if the 
hydrodynamic interactions between the particles are neglected. The divergence at zero separation is 
the signature of the random attractor asymptoted by the particles' trajectories at large times. How- 
ever the hydrodynamic interactions produce a repulsion between the particles that is non-negligible 
at small separations. We introduce equations governing the repulsion and show it smoothens the 
singular attractor near the particles where the pair correlation function saturates. The effect is most 
essential at the Stokes number of order one, where the correlations decrease by a factor of a few. 
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The effect of inertia on the motion of particles driven 
by the stationary turbulent flow has enjoyed much atten- 
tion of the researchers recently p]4l2j. The progress in 
the understanding of the behavior of tracer particles [l3[ 
made it natural to try to understand the impact of inertia 
that real particles have always. It was recognized quite 
early that inertia brings new qualitative effects Q and 
even small inertia can have a profound effect on the spa- 
tial distribution of particles in the flow . However small 
the inertia is, the particles' motion in real space has a ran- 
dom attractor: at large times the trajectories asymptote 
a multi-fractal set that evolves in time keeping its sta- 
tistical properties constant [4, 5]. The full set of fractal 
dimensions of the attractor was found recently [H 03] 
for real turbulence. The derivation involves no modeling 
and considers the statistics of turbulence as given but 
unknown. The statistics of density is log-normal and de- 
termined completely by the pair-correlation function [3] . 
The latter obeys a power-law with a negative exponent 
at small separations pj-[l2j]. The divergence at zero sep- 
aration is the signature of the singular spatial structure 
of the attractor, while the exponent gives the correlation 
codimension of the attractor. These results are derived 
neglecting the hydrodynamic interactions, while one of 
the main applications is to collisions where the particles 
come near each other and the interactions become im- 
portant. An example is the formation of rain in clouds 
where droplets can grow thanks to the collisions with 
other droplets, see 0, Ell and references therein. The 
rate of the collisions is proportional to the probability 
that turbulent fluctuations bring the particles together 
and it is measured by the pair-correlation function of the 
density at the separation equal to the diameter of the 
particles. The hydrodynamic interactions produce a re- 
pulsion that can be strong enough to prevent the collision 
completely. An example is the fall of a heavier particle 
onto the lighter one in the still air under the action of 
gravity. The repulsive force grows as the inverse of the 
distance between the particles' surfaces, and in the frame 
of hydrodynamics the particles never collide [l5l-[l9j . 



In this Letter we derive equations that include the hy- 
drodynamic interactions and describe the effect of the 
interactions on the pair-correlation function of particles' 
density. The analysis is performed at small inertia, as 
measured by the dimensionless Stokes number St. The 
depletion of the power-law in the pair-correlation func- 
tion due to hydrodynamic interactions is found. The 
results are extrapolated to assess the impact of hydro- 
dynamic interactions at St ~ 1, where the preferential 
concentration is maximal 0. 

The hydrodynamic interactions between two particles 
do not admit exact analytic expressions and a simpli- 
fication is necessary to perform the theoretical analysis 
[l5l - [20| | . For heavy particles considered here the hydro- 
dynamic interactions are determined by the perturbation 
flow caused by the inertial slip of the particles with re- 
spect to the surrounding flow. This flow obeys linear 
equations and the method of superposition is available 
to derive the forces acting on the particles as sums of 
special cases. The hydrodynamic interactions are impor- 
tant mainly when two particle move toward each other 
along the line of centers with equal velocities dU . Say, 
when the distance between the surfaces of the particles 
measured in the units of their radius is 0.01 the differ- 
ence between the forces in this case and all other cases 
is about two orders of magnitude [H, Gj| ■ To capture 
the main effect we account only for the hydrodynamic 
interactions caused by the motion of the particles toward 
each other. This is reasonable as this is the component 
of the motion describing the particles' approach. 

We assume that without particles the forces stirring 
the fluid produce the turbulent flow u(t,x) and the pres- 
sure po(t, x) that solve the Navier-Stokes (NS) equations 
[2ll |. The flow in the presence of particles v(t,x) obeys 



d t v + v ■ Vv = -Vp + vV 2 v, V 



0, 



(1) 



where the particles are accounted via the no-slip bound- 
ary condition (b. c.) v(\x — Xi(t)\ = d/2) = Vi(t) and 
far from the particles the flow v and the pressure p must 
match u(t,x) and po(t,x). Here and Vi are the par- 
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tides' coordinates and velocities, respectively. The par- 
ticles' diameter d is assumed to be much smaller that 
the Kolmogorov scale of turbulence rj determined by the 
kinematic viscosity v, see [2lJ. The flow around particles 
separated from the rest by distances much larger than d 
is the sum of u and the Stokes flow caused by the rela- 
tive motion of the particle through the flow at the speed 
Vi — u(t,Xi). For heavy particles the latter is damped 
via the linear friction force TVi = u[t, Xi(t)] — Vi where r 
is the Stokes' time and the Reynolds number i?e„associ- 
ated with the perturbation flow is assumed small [l| . The 
efficiency of relaxation is measured by the Stokes number 
St = At, where A ~ v/rj 2 is the characteristic value of 
|Vii| and A -1 is the minimal time-scale of variations of 
u. At St <C 1 the particles follow the flow closely and 
one finds Vi = u — r[dtu + (u ■ V)u], where the RHS 
is evaluated at a^. The estimate for Re p based on the 
latter expression validates Re p <C 1 at St <C 1. Thus at 
St -C 1 one can introduce the particles' velocity field 



x i = V[t,x i (t)], V = u-r[d t +«• V]u. (2) 

While V — ti is relatively small, it brings a qualitatively 
new effect: V • V = —r\/jUiViUj ^ 0, though V • u = 0. 
The continuity equation <9 t n + V • (riV) = on the par- 
ticles' density n has no constant solution at r > 0. The 
particles trajectories approach a multi-fractal set in space 
and the steady-state fluctuations of the density obey 

(n(O)n(r)) = (tj/r)"", r « V , n = 2 |^ Ai/A rf | , (3) 

where the angular brackets stand for spatial averaging 
[3l-fl2l|. (n) = 1 and Xi are the Lyapunov exponents 22 1 
of the particles' flow V. Both an d A^ are neg- 

ative and their leading order behavior at 5t < 1 is 
| J2 A, | oc St 2 and A d « \ d (St = 0) = A$, where AJ£ ~ A 
is the d— th Lyapunov exponent of u. The latter de- 
termines the exponential diver genc e of fluid particles' in 
turbulence back in time [IGJJll- Importantly for the 
following, d t u contribution into v — u in Eq. (J3J has zero 
divergence and can be neglectedj_as can be seen explic- 
itly from the complete solution [12]. Eq. ([3|) is derived 
analytically Q at St <C 1. As St increases, the power- 
law holds for St < St cr i t i ca i where n first increases and 
then decreases to become zero at St = St cr itioai- At 
St > St cr i t i ca i the power-law breaks down @. The pref- 
erential concentration is maximal [9] at St ~ 1, but sig- 
nificant @ already at St <C 1. This is because St — » 
is a singular limit: linv^o lim5 t _ >0 (n(0)n(r)) = 1, while 
limst_>o ]ha. r -to{n(0)n(r)) — oo, so at a sufficiently small 
scale the inertia is always important. 

The hydrodynamic interactions demand a different 
framework: the particle's velocity is determined not only 
by its position, but also by the position of the particle 
with which it interacts, so the particles' flow is no longer 
defined. To find (n(O)n(r)) we use that it is proportional 
to the PDF P(R) of the distance R between two particles 
in the steady state. The existence of the latter demands 
a cutoff R cu t at large R, as on average the two particles 



disperse. The cutoff does not enter the final result after 
the proportionality constant is fixed by the demand that 
(n(O)n(r)) tends to (n) 2 = 1 at large r. Say if the cutoff 
is due to a finite volume, then the particles spend most 
of the time at separations comparable with the volume's 
size. Occasionally they are brought close by the random 
flow. The frequency of these events is measured by P(R). 

We first re-derive (n(O)n(J?)) from P(R) at R ^> d 
where the hydrodynamic interactions are negligible and 
one can use Eq. ([2]). At St <C 1 or /i <C I there are no 
fluctuations of n at R ~ rj, see Eq. ([3j), so we assume 
R<r). Then the distance obeys R=V[t, x(t) + R(t)} - 
V[t,x(t)] m crR. Here = \7jVi{t 7 x = x(t)] can be 
expressed via = S7jUi[t,x = x(t)] as a — s — t[s + s ], 
see Eq. ([2]). By incompressibility trs = 0, so tro = —trs 2 . 
One has R{t) = W{t)R(0), where W = aW, so 

P(R) = (5 [R(t) -R}) = (6 (R(0) -W~ l R) / dct W) . (4) 



Introducing £ 



= rtrs 2 one may write detW(t) — 
exp(— ((t')dt'). From the definition |J^Aj| = 
-lim^oolndet W(t)/t we have (Q = \ J2M> see HI- 
The "Fluctuation-Dissipation Theorem" ("FDT") holds 
(C) = (1/2) f{{((Q)((t)))dt, where the double angular 
brackets designate dispersion [23| . At t much larger than 
the correlation time A -1 of a, see [2lJ, one can neglect 
in Eq. (U) the contribution of the time interval of order 
t c near t = and perform independent averaging over 
W(t) and the velocity field at t < 0. The latter av- 
eraging averages -R(O) producing the steady state condi- 
tion p(r) = <jex P ( r 4 c(£')^') p{w- i r) x j : cf. mm. 

The solution (independently of R cu t) is P(R) oc R~ a , 



exp 



((t')dt') W-\t)R 



1, R=R/R. (5) 



At St < 1 one has « « 1 giving \W~ 1 {t)R\~ a w 
cxp[— a|A^|t]. Furthermore at St -C 1 one can use Gaus- 
sian approximation for (exp ( f£ ^(t')dt'^) which using 

the FDT gives (exp ( J* ((t')dt'^)) w exp (2t\ J2 A»|). The 

substitution in Eq. reproduces a = fi in Eq. ©. 

The observation that the divergenceless part of the in- 
ertial correction V — u can be omitted from V without 
affecting the results to the leading order in St implies 



R = a'R = sR - (R/D, ct- = 



SijC/D, (6) 



where D = 2, 3 is the dimension, should give the same a 
(the part of a describing the divergence of V is (Sij/D). 
This can be verified substituting W — W in Eq. ((4]) 
where W = a'W . Equation (|6|) says that outside the 
vortices, at trs 2 > 0, the particles are attracted due to in- 
ertia. Inertia produces a spherically symmetric contrac- 
tion or expansion of all R, impossible for incompressible 
flow, while the rest of its effects are negligible. 

Using Eq. §5§ with W = W one also sees that the 
statistics of s in Eq. (O can be considered as the one of 
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V jUi in the fluid particle's frame and not the particle's 
frame, while s can be assumed independent of £. The 
values of | W — that determine the average in Eq. JSJ) 
are exp(|A^|£), so P(R) is determined by the events for 
which R decreases exponentially from the most proba- 
ble, large, values of R at the rate |A$|. This is the same 
for fluid particles as |A^| determines the rate of exponen- 
tial divergence of trajectories backward in time. Inertia 
changes the volume of points that contract at this rate. 

To account for hydrodynamic interactions, we consider 
two particles separated by R ~ d, well-separated from the 
rest. We look for the flow around the particles in the form 
v(x) = u(x) + U[x — x(t)] &ndp(x) = p (x) + P[x — x(t)]. 
Substitution in Eq. (UJ and the use of Re p <C 1 give, 

VP = vV 2 U 1 U(x = d/2) = x-u[t,x(t)], (7) 
U(\x -R\ = d/2) = x + R -u[t, x(t) + R(t)} ,(8) 

with U, P vanishing at large x. This is the problem of 
two particles moving at speeds J\ = x — u [t, x(t)] and 
J2 = x + R — u[t,x(t) + R(t)} in the fluid at rest at 
infinity (lij . The motion can be written as (i = 1,2) 

J t = (-iyj r R/2 + J t - (-l)V r jR/2, J=J 2 Ji-(9) 

where J r = J R. We consider the forces Fi on the 
particles as the superposition of the forces produced by 
the first term and the rest. At J r < the former de- 
scribes particles' motion toward each other along the 
line of centers at the same speed J r /2. In this case 
the particles experience equal resistance F(R)\J r \/2r (we 
set the mass to unity). At R d Stokes' law holds, 
F(R) « 1. For R ~ d the flow is not the superposi- 
tion of two Stokes flows and F(R) ~ 1/(1 — d/R), see 
[l5j . This behavior makes the collisions in the still air 
impossible pj|. We use F(R) = 1/(1 - d/R) for quali- 
tative estimates. The rest of the terms in Eq. © do not 
change the distance between the particles and the de- 
viation from Stokes' law for them is much weaker [l5j . 
Studying the impact of hydrodynamic interactions on 
the particles' approach to each other, we account for 
the main effect and neglect the latter deviations, so that 
tF, = -Ji - {-l) i e (- J r ) J r [F(R) - 1] R/2, where 9(x) 
is the step function. Subtracting the equations on x = F\ 
and x + R = F2 and considering R -C r\ we find 

tw = sR-w+[1-F (R)] (w r -RsR\ 9 RsR-w r R, 

where w = R. At St — one finds the equation for fluid 
particles w = sR: there are no interactions for inertia- 
less particles that cause no friction. Solving the equation 
to the first order in St and using a — s — t[s + s 2 ] we find 

R = aR + zO (z) [1 - F- 1 (R)] R, z = tR[s + s 2 ]R. 

To understand the effect of the hydrodynamic interaction 
force, which is radial, we multiply with R. Introducing 
p = \n(R/d) and £ = RsR we have p = £ - z/F(R) 



at z > 0. Compared with p = ^ — z holding without 
interactions, the hydrodynamic interactions deplete the 
inertial correction to the radial component of the relative 
velocity by the factor of F(R), when the correction de- 
scribes approaching particles. Thus the interactions are 
expected to eliminate the inertial enhancement of P{R) 
at R ~ d, so P(R) and (n(0)n(R)) saturate at R ~ d. 

The effective description explained previously account- 
ing only for the trace of s + s 2 gives (F^ 1 = 1 — d/R) 



R = sR~ [l-d/R](R/D, 



(10) 



for £ > and Eq. © otherwise. Equation (fTU)) trans- 
forms to Eq. ([6]) at R ^> d and thus reproduces the results 
at R 3> d, while capturing the main effect of the hydro- 
dynamic interactions at R ~ d: the cancelation of the 
approaching inertial component of the relative velocity. 
The particles can collide due to sR term but the inertial 
enhancement of the collision velocity vanishes. Equation 
(1101) should be supplied with the boundary conditions (b. 
c.) at R = d, describing the physics of the collision, e. 
g. absorbing b. c. for coalescence. Consider the events 
that form the inertial enhancement of P(R) at R ~ d 
and thus have trs 2 > 0. As one tracks R(t) back in time 
it grows exponentially at the rate \Xd\- The interactions 
become negligible at t ~ — | | 1 where R(t) ^> d. Since 
trs 2 varies over the same time-scale | | 1 , then the be- 
havior of P(R) at R ~ d should be captured by assuming 
trs 2 > during all the time-interval the interactions are 
relevant. Thus to understand the behavior of P(R) at 
R ~ d we use Eq. (fTU]) for any sign of C- We find 

p = £-(l-exp[-p])C/A R = sR-£R. (11) 

The dynamics of the orientation decouples and £in the 
equation on p can be considered as given, cf. [251 ] . As 
mentioned, s (and thus £) and £ in can be considered 
independent where the statistics of Sjj can be considered 
as the one of V ' jUi in the fluid particle's frame. Then £ is 
the same noise that governs the exponential separation of 
fluid particles, so e. g. (£) is the first L yap unov exponent 
of the motion of fluid particles, cf. [11 H3, H|. The 
Kraichnan model [l3| prescribes £ = A" + £' where £' 
is a white noise (£'(ii)£'(* 2 )) = 2A?<5(* 2 - h)/D. This 
and C = I I + C define the Kraichnan model for the 
problem, where (/ is a white noise which amplitude is 
fixed by the "FDT" as (C'(*i)C'(*2)) = 21X^1^2 ~ 
ti). One can write the Fokker-Planck equation on the 
PDF P(p) of p which gives that in the steady state the 
probability current must be constant. The constant is 
zero for reflecting b. c. at p = 0. This gives the condition 

P' _ 2D 2 - pD + p(D - 1) exp[-p] + ^exp[-2p] 
P ~ 2D + p{l-exp[-p}) 2 

where p = 2| ^ A$|/|Ad|, see Eq. ©, and we used [l3[ 
that in the Kraichnan model Ai = |Ad|. Expanding the 
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FIG. 1: Shown is the plot of log (P(p) exp[— 2p\) versus p. 
The power-law starts saturating at p ~ 2.1 for St — 0.4 and 
p « 2.4 for St = 0.2. 



RHS to the first order in p <C 1 and integrating, 

P = Nexp^pD - pp + p(3D - 1) (1 - exp[-p]) /(2D) 
-/x(P>-l)(l-exp[-2p])/(4£)) X ), (12) 



where iV is the normalization factor. At p — we recover 
P{p) oc exp[pD] equivalent to P(R) = const. The hydro- 
dynamic interactions produce the last two terms in the 
exponent. At exp[— p] < 1 or r > d these saturate at a 
constant and (n(O)n(r)) cx P[ln(r /d)]{d/r) D reproduces 
Eq. ([3]). At p <C 1 the linear term in the Taylor ex- 
pansion of the hydrodynamic interactions' terms cancels 
-HP and (n(O)ra(r)) « A>exp[p(P + 1) ln 2 (P/,i)/(4P>)] 
depends on r weakly. Between p < 1 and exp[— p] -C 1 
there is a smooth transition between the two behaviors. 



These modifications are less important at p <C 1 where 
(■q/lQd)^ « {r]/d)^ and {rj/r)^ near the particles is almost 
the same as far from them. In contrast, at St ~ 1 where 
Eq. (J3|) holds with /x ~ 1, the effect is important and 
(n(0)n(d)} is lowered by the factor of a few, see Eq. (fT2|) . 

The above features are expected to be generic. We 
verify this numerically for a model where a has a corre- 
lation time of order of A -1 , mimicking turbulence. We 
consider the D — 2 random renewal model [25| where the 
traceless, zero mean matrix s a p(t) is a piecewise constant 
process, s(t) — s k at kAt < t < (k + l)At generated by 

s k ap = C(At)- 1 [ e ^ 3 V* + /*/^/V2-^(/ fc ) a /2V2], 

where e Q ^ 7 is the antisymmetric symbol, V k and / * are 
independent standard Gaussian variables. At C — > 
one gets the Kraichnan model (25|. We consider C ~ 1. 
We model ( as | Yl Ai I + C where the zero- mean Gaussian 
noise £' is renewed each At and obeys | $^ Aj| = (C 2 ) At/2 
to ensure the "FDT". The results of 4 x 10 6 renewals 
for C = 2, St = 0.4 and C = 1, St = 0.2, where 
St 2 = \ J2Xi\At are shown in Fig. [TJ When p reached 
p = it was reset at p cut = 1800 (underestimating P(p) 
near p = 0). At p = p maa; reflecting b. c. were used. The 
results agree qualitatively with the Kraichnan model. 

We introduced the equations describing the two- 
particle dynamics of inertial particles at small St with 
the account of hydrodynamic interactions. The interac- 
tions eliminate the inertial effects in the particles' vicinity 
and the divergent power-law in the correlation function 
of the density saturates. The density is smooth at scales 
r ~ d, in contrast to the singular density in the single- 
particle approximation. Our equations are useful even 
without the interactions, for problems such as collisions 
of water droplets in clouds mediated by turbulence. 
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